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1. The problem 

In a cosmo logical context, where there is a preferred four- velocity field u a , the electric 
and magnetic parts of the Weyl tensor, E ab and H ab , represent covariantly the locally 
free gravitational field, in the sense that they are not pointwise determined by the matter 
fields (see e.g. 0-0). They do not arise explicitly in the Einstein equations. Although 
they are not really independent of the matter fields, being constrained by integrability 
conditions in the form of the Bianchi identities, it is nevertheless useful to think of them 
as being locally free. 

The Bianchi identities involve the time derivatives and the spatial divergences and 
curlsf of E ab and H ab (see the Appendix). The locally free parts of the curvature are not 
pointwise determined by the matter fields via Einstein's equations. Similarly, the locally 
free parts of the derivative of the curvature are those parts not pointwise determined via 
the Bianchi identities by the matter fields and their derivatives and by E ab , H ab . The aim 
of this paper is to determine these locally free parts and their relation to gravitational 
waves in cosmology. 

We find via a complete covariant decomposition that the parts of V ' c E ab and V c H ab 
which are missing from the Bianchi identities are exactly the projected totally symmetric 
and trace-free parts V{ c E ab } and V( c #a&)4 We call these parts the distortions of E ab 
and H ab . The distortions represent the part of the derivative of the curvature which 
is uncoupled from the Bianchi identities and does not locally interact with the matter 
fields. The divergences are determined pointwise by the matter terms and the locally 
free field. The curls (together with the matter terms and the locally free field) determine 
the time derivatives of E ab and H ab pointwise. Therefore the distortions and the curls 
provide a covariant characterisation of the locally free part of the space-time gradient of 
the gravitational field. 

We show that in the case of gravitational waves (covariantly characterised), the 
Laplacian is determined by the distortion. Hence not only the curls, but also the 
distortions of E ab and H ab must be non-zero if there are gravitational waves. This 
result represents a refinement of the covariant characterisation and understanding of 
gravitational waves in cosmology (see e.g. || - W% )- 

In Section 2 we give the complete covariant decompositions of the derivatives of the 
Weyl tensor and fluid kinematic quantities. The derivations of these decompositions, 
involving a combination of covariant tensor methods and Young diagram methods, are 
presented in Section 3. In Section 4 we show how the distortions and curls determine the 

f The terms divergence and curl as applied to rank-2 tensors are defined by equation ([!]). 

\ Angle brackets (• • •) enclosing indices represent their totally symmetric trace-free and spatially 

projected part. Round and square brackets denote, respectively, symmetrisation and anti- 

symmetrisation. 
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covariant gravitational wave equations. Further implications of our results are discussed 
in Section 5. The Appendix gives the Ricci and Bianchi identities in 1 + 3 covariant 
form, together with the required identities for second order derivatives. 

2. The basic results 

The locally free field E a b, H a b and its characterisation are covariant if a covariantly 
defined four-velocity exists. For a single perfect fluid, there is a unique covariant four- 
velocity, whereas for imperfect fluids and multi-fluids there is in general more than one 
covariantly defined four- velocity. For simplicity, we restrict ourselves to the perfect fluid 
case. The results are readily extended to the imperfect fluid or multi-fluid case, given 
any covariant choice of four-velocity. 

For a gravitational field with perfect fluid source, the basic covariant variables are: 
the fluid scalars 6 (expansion) , p (energy density) , p (pressure) ; the fluid spatial vectors 
ii a (4-acceleration), uj a (vorticity); the spatial trace-free symmetric tensors a a b (fluid 
shear), E a b, H a b] and the projection tensor h a b = g a b + u a Ub, which projects orthogonal 
to the fluid 4-velocity vector u a . The covariant spatial derivative is defined by 

D c S a - b ... = h c r h%---h b q ---V r S p - q .... 

It follows that D c h a b = 0. The covariant projected permutation tensor is e a b c = rj abc dU d , 
satisfying D a e 6crf = and e a b c e de ^ = 3\h[ a d hb e h c ^ . The covariant spatial divergence and 
curl of rank-2 tensors are defined by |J 

(div S) a = V b S ab , (curl S) ab = e cd{a D c S b) d , (1) 

generalising the expressions for vectors: 

div V = D b V b , (curl V) tt = s abc D b V c . 

For brevity, we shall henceforth omit the brackets around 'curl' unless this leads to 
ambiguity. 

The aim now is, through the covariant 1 + 3 splittings of vectors and rank-2 and 
rank-3 tensors, to find the complete covariant decomposition of the derivatives of scalars, 
vectors, and rank-2 symmetric tensors. Any space-time vector S a may be covariantly 
split as: 

S a = -S b u h u a + h a c S c . 
As a special case, the derivatives of the fluid scalars have the decomposition 

v a e = -e Ua + D a e , 

and similarly for p and p. 
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Any spatial tensor S ab = h a c h b d S cd may be covariantly decomposed as: 

Sab = S{ab) + \h c dS cd h a b + e a b c S c , (2) 

where S( a b) = h^ a c h b ) d S c d — \h c dS cd h a b is the projected symmetric trace-free part, and 
the spatial vector S a is equivalent to the skew part: 

S[ab] = ^abcS "vv - S a — \s a bcS • (3) 

If S a b = S( a b), then (curlS 1 )^ = (curl S") ( a h) . If S a b = S[ a b\, then we find the identities 

curl S ab = V {a S b ) - p c S c h ab , D b S ab = curl S a . (4) 

The decomposition of the derivative of the fluid 4-velocity [|IJ is in our notation 

V 6 w a = -u a u b + D 6 w a , (5) 
Db«a = a ab + l<S>h ab + uj c e ab c , (6) 

<?ab = D(aMfe) , @ = D°M a , UJ a = -^Cm\u a , (7) 

where (6) is a particular case of @. For the vorticity, using @ - (|J) we find 

VfMj = -U c U C U a U b - LO(a)U h - U a {ctbc^ + + Dfe^a , (8) 

D b uj a = D (a ^ 6) + lB c io c h ab - lcurluj c €ab c , (9) 

where S( a ) = h a b Sb, and similarly for the 4-acceleration ii a . 

The electric and magnetic parts of the Weyl tensor, E ab = C ac bdU c u d and H ab = 
^£acdC cd b e u e , satisfy E a b = E^ a b) and H a b = H( a b)- Their derivatives have the following 
covariant decomposition: 

V c E ab = - E {a b)U c - 2u( a E b ) d u d u c 

+ 2u {a [E b) d (a dc + u dc ) + \QE b)c } + D c E ab , (10) 

B c E ab = E^ + ~D d E d{a h b )c - fcurl E d{a e b)c d , (11) 

where S( a b) = h( a c hb) d S cd — \h c dS cd h a b, and where E ca b is completely symmetric and 
trace-free: 



E C ab = E^ab) = ^>{ C E a b) • (12) 

du £ 
{a lib) 



Note that the projected time derivative in (10) satisfies h a d h b e E de = h/ a d h b \ e E de . This 



decomposition is derived in the next section. 



A similar decomposition applies to V c H ab . The generalisation (11) of (6) suggests 
that we call E ca b and H cab the distortions of E a b and H a b, by analogy with the shear o a b- 
These represent the divergence-free and curl-free spatial variation of the Weyl tensor. f 

f Note that the term distortion in the case of a rank-2 tensor is not intended to imply the geometrical 
significance that it has for a vector Q]. 
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The decomposition (10) - (|T^) applies to the derivative of any spatial, symmetric and 
trace-free rank-2 tensor; thus we have for example 

D c 0a& = o cah + ^D d a d{a h b)c - |curl a d{a e b) d , (13) 

in the case of the shear a ab . 

Using the notation of ||, the covariant splitting of the Ricci and Bianchi identities 
given in |1] is considerably simplified. For convenience, these equations are given in 
the appendix, in the case of a perfect fluid source for the gravitational field. All the 
covariant time and spatial derivatives of the scalar and vector variables - ©, p, p, u a , 
uj a and u a - appear in the covariant Ricci and Bianchi equations. This includes the 
distortions D{ b u a ) = a ab , D/ b a> a ) and D(pii a \. But for the tensor variables E ab , H ab 
and a a b, the distortions do not appear in any of these equations. In the case of the 
shear, the divergence and curl terms occur in the Ricci constraint equations (A4) and 
(A6) respectively, but the distortion a cab does not occur in any of the Ricci or Bianchi 
identities. However, the evolution of a cab is governed by the distortion of E ab , as may be 
seen by taking the symmetric and trace-free part of the spatial derivative of the shear 
propagation equation (A3). 

The distortion arises in the covariant Laplacian, which is central to the wave equation 
(see Section 4). Consider the distortion of say H ab : 

H C ab = ^( c H ab ) — | h(abD H c )d . 

Taking the divergence, we get 

^> c H cab = ^D 2 H ab — ^D( a D c if 6 ) c + |D 6 D^ a if 6 ) c , 

where D 2 = D a D a is the covariant Laplacian. The last term on the right may be 
converted into a divergence term plus curvature correction terms, using the commutation 
identity (A13): 

B 2 H ab = 3B c H cab - 2 (p - f 6 2 ) H ab - p {a D c H b)c + M[H] ab , (14) 

where 

M[H} ah = 4u c uj c H ab + 2Qa c {a H b)c - 6u c u {a H b)c 

+ 2a c ( a o b ) d H cd — 2a cd H cd a ab + 2Qu c e cd i^ a H b ) d 
— 2o cd o c ( a H b ) d — <oE c ( a H b ) c + 2uj c e cde (j d ( a H b } e 
+ 2u c e cd ( a a b ) e H de + 2uj c e cd a H^ bd ) + 2u c e cd b H( ad ) . (15) 

The new identity (0), (15) shows explicitly how the covariant Laplacian is determined 
by the divergence of the distortion, the gradient of the divergence, the time derivative, 
and algebraic terms. The homogeneous and isotropic Friedman-Robertson- Walker 
(FRW) universe is covariantly characterised by the vanishing of u a , u a and a abl which 



6 



then forces E a b, H a b and the spatial gradients of p, p and to vanish. In an almost 
FRW universe (which is considered in Section 4), these covariant quantities are small 
||, and the term J\f[H] a b is non-linear and may be neglected. 



3. The decompositions 



The derivations of the covariant decompositions (5) - (9) for the vectors u a and uj a 
(and the corresponding expressions for u a ) involve a straightforward application of the 
1 + 3 splitting of rank-2 tensors together with @ - @). For the tensor derivatives, we 
need the 1 + 3 splitting of rank-3 tensors, the generalisation to rank-3 tensors of the 
decomposition (0), and the properties of the covariant tensor curl and divergence. 
A: The space-time splitting. Any rank-3 space-time tensor has the covariant 1 + 3 
decomposition 

Sabc = au a u b u c + V a u b u c + W b u c u a + Z c u a u b 
+ A a bUc + B bc u a + C ca u b + F abc , 



where V a , ■ ■ ■ , F abc are spatial tensors. Since E ab 
u a u b E, 



E, 



(ab) 



and E ab u b = 0, we have 



ab:c 



and so 

E a b;c = 2V {a u b) u c + A ab u c + 2u {a B b)c + 

Fabc 



where A ab = A^ ab ) and F abc = D c E ab . Contracting with u c and u b and using (6) gives 
A a b + V a u b + V b u a = -E^, 

B ac + V a u c = (\®h b c + a b c + u b c - u b u c ) E ab , 

and (10) follows. 

B: The spatial splitting. Any rank-3 tensor in n dimensions may be decomposed as 
M abc = M( afec ) + S abc + T abc + M[ abc ] , (16) 

where 

Sabc = \M[ab\ c + |M[ | c |b] = — S ba c > 

Tabc = lM a [ bc ] + §M[ fe | a | c ] = — T ac b ■ (17) 

In terms of Young diagrams [12|], this is represented by the irreducible decomposition 
O O |~c~| = ^ 



b 


c 


a 





b 


a 


c 





so the 4 terms on the right of (16) are linearly independent. For n 
10 + 8 + 8 + 1 = 27 independent components. (In general, they have 

±n(n + l)(n + 2) + \n{n + l)(n - 1) + \n(n + l)(n - 1) 



3, they have 



+ ±n(n - l)(n- 2) 



n 
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free components.) Since the terms in (|T^) are linearly independent, we can split off their 
trace-free parts independently. Because of the symmetry, each has (at most) one trace 
term. Using the appropriate metric hob- 

M {abc) = M {abc) + — — M {a h bc) , (18) 
n + 2 

2 

Sabc = Sabc H -S\ a h b ] c , (19) 

n — 1 

2 

Tabc = T a b c H -T[bh c ] a , (20) 

n — 1 

where M a = M( a bc)h bc , S a = S abc h bc , T a = T abc h ac , and S abc , T a bc are trace-free. If 
M a b c = M(ab) c , then symmetrising fllBD on ab gives 

M abc = M(abc) + ^(ab)c , (21) 

and (17) becomes 

Tabc = f M a [6c] . (22) 

Notice that, when M a fe c = Mb ac , anti-symmetrising (|2lD on be gives, using (|22|), 
|T abc = T (afe ) c - T (ac ) b , so T (ab)c is equivalent to T afec . As before, we can split off the 
trace-free parts of (^T|). If M a bc = M[ a b]c, then T| a b] c = and (p!6|) becomes 

Mabc = S a bc + M[ abc ] ■ (23) 

Specialising now to the case of spatial rank-3 tensors, so that n = 3, simplifications 
result from (^|): 

M[ a bc] = MSabc , Sabc = Sdc^ab 1 , ?afcc = TdaSbc , (24) 

where S cd = \S a bc£ ab d = S{cd) and similarly for f cd . If M abc = M (o6)c , then by (18), (20) 
- (|22]) and (^4|), we have 

Mibc = M( ab )c = M^bc) + |M( a /i fec ) + \T d ( a £b)c d + \T[bh c ] a + T[ a h c ]b , (25) 

where M a = \Mj - \M d d a and T a = \M ad d - \M d d a . If in addition M a \ = 0, then 
M a = T a and (p5|) becomes 

M a6c = M ( ab)c = M {ab c) + ^M {a h b )c + f d(a Eb)c d • (26) 

If M abc = M [ab] c, then (0) - @§ give 

Mibc = ^[afe] c = S[ a hb] c + S dc e a b d + Me a bc , (27) 
where S" a = M afc fe . 

Equations (|25|) - (|27|) give the complete covariant decomposition of any rank-3 
spatial tensor. An example of (|27D is provided by the decomposition of the commutators 
Yab of an orthonormal triad fI5|| . The decomposition ( |26|) applies to D c E a b, D c H a b and 
~D c (J a b- In the case of B c E a b, contracting (p6|) with /i 6c gives D b i? afe = |M a , while 
contraction with e cb e and symmetrisation on ea gives curli? ea = — |T ea . Thus (11) is 
established. The corresponding relations for H a b and o~ a b follow immediately. 
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4. Gravitational waves 



The identification of the distortions as the parts of the Weyl derivative that are 
uncoupled from the Bianchi identities, suggests that they play an important role in 
the existence of gravitational waves. We show here that this is indeed the case. 

Gravitational waves in cosmology are gauge-invariant transverse tensor perturba- 
tions of the homogeneous and isotropic FRW spacetime. The covariant characterisation 
of gravitational waves in terms of the Weyl tensor was introduced by Hawking || and 
developed in the covariant and gauge-invariant perturbation theory of Ellis and Bruni 
||. In this approach, the fundamental equations are the covariant propagation and 
constraint equations (Al) - (A12) given in the Appendix, with the right hand sides set 
to zero (i.e. linearised). Gravitational waves are described by E ab and H ab , governed by 
the linearised form of the Bianchi identities (A9) - (A12). In the absence of scalar and 
vector modes, the gradients of p, p and 6 are zero, together with u a and u a . It follows 
from the linearised version of equations (A4), (All) and (A12) that a ab , E ab and H ab 
are divergence-free. These tensors satisfy wave equations ||, [|7| provided their curls 
are non-zero (strictly speaking, the curl of the curl must be non-zero). For example, 
taking the time derivative of the linearisation of equation (A10), using the linearisation 
of equations (Al), (A6), (A9) and the identities ( |A14| ) and ( |A15| ), we get 

n2 H ab = — H ab + B 2 H ab 

= leH ab +(l& 2 -2p)H abj (28) 



in agreement with J7[ and p4| . 

The derivation of (28) reflects the known result J7| that a necessary covariant 
condition for gravitational waves in cosmology is 

D b E ab = + curl E ab , D b H ab = + curl H ab . (29) 

Now the identity (^|) derived in Section 2 for the covariant Laplacian allows us to refine 
this result. In the linearised case, the term M[H] ab given by (15) vanishes, and we have 



D 2 H ab = 3D c H cab - 2 (p - |6 2 ) H ab , (30) 

since the divergence is zero for gravitational waves. It follows that if the distortion of 
H ab vanishes (strictly, if the divergence of the distortion vanishes), then H ab satisfies a 
covariant Helmholtz equation 

D 2 H ab = -(^jH ab , (31) 

where we have used the background Friedmann equation to evaluate the coefficient in 
(PD|), with a the scale factor and /C = 0, ±1 the spatial curvature of the background. 
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For K, = 0, the Laplacian of H ab vanishes and there can be no wave propagation. For 
JC 7^ 0, H ab is proportional to a tensor eigenfunction of the covariant Laplacian. These 
covariant eigenf unctions satisfy || 

1,2 

n 2 n( fc ) — r>( fc ) r>W — r> b r>( k "> — n — n( fe ) 

where /c determines the comoving wave number. For JC — — 1, (|31"D shows that there is 
no real value of k (and thus no smooth eigenfunction), which forces H ab = 0, so that 
there are no gravitational waves. For JC = +1, (^1|) implies k = \/6, so that there is at 
most one wavelength admitted. However, this rules out general gravitational waves. A 
similar result follows if we start with zero distortion of E ab . Thus we have the extension 
of ©: 

a necessary covariant condition for the existence of gravitational waves in cosmology is 
D b H ab = 0, cur\H ab ^0^D {c H ab) , 

D b E ab = , curl E ab ^ ^ B {c E ab) . (32) 
5. Conclusion 

We have found the complete covariant decomposition of the derivatives of the electric 
and magnetic parts of the Weyl tensor. This extends the covariant characterisation of 
the locally free gravitational field to its space-time gradient. The derivatives decompose 
into time derivatives and spatial curls, divergences and distortions (totally symmetric 
trace- free parts). The covariant Ricci and Bianchi equations given in the Appendix 
involve all these parts of V c E ab and V c H ab except for the distortions E cab = D^Ea^ 
and H cab = D( c H ab ), which do not occur in the equations. They may be thought of 
as the uncoupled parts of the derivatives, containing at most 2 x 7 = 14 independent 
components out of the maximal 40 independent components of V ' e C a bcd- Together with 
the curls, the distortions are crucial for the existence of gravitational waves. 

At any spacetime event V, specification of C abc d and V ' e C a bcd determines a first-order 
Taylor approximation to C abc d near V . The divergences are determined algebraically at 
V via the Bianchi constraints (All) and (A12). The curls algebraically determine the 
time derivatives via the Bianchi propagation equations (A9) and (A10), so determining 
the future evolution of the system. Together with the distortions at V, they are free 
data at that point, in that they are not determined by the matter there (i.e. by p, p, u a 
and their first derivatives). Rather their values there are determined non-locally, their 
time evolutions along the fluid flow-line through V being determined by the free data 
at earlier times (through spatial derivatives of the Bianchi propagation equations (A9) 
and (A10) for E ab and H ab ), and their spatial derivatives being restricted by spatial 
derivatives of the Bianchi constraint equations (All) and (A12). 
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It is clear from this that the distortions at V are the least coupled parts of the 
derivative of the Weyl tensor. Together with the curls, they are the free data that 
can be given at a point. Since the curls contain at most 2 x 5 = 10 independent 
components, there are at V in general 14 + 10 = 24 unconstrained components of 



V 'eCabcd- These tensor properties at a point are reminiscent of Penrose's result [FJ| that 
the totally symmetrised derivatives of the Weyl spinor are freely specifiable at a point 
in the vacuum case (see []nj for the non-vacuum generalisation, and ]TB| for further 
discussion and application, of Penrose's result). 

We have shown how the distortions of H ab and E ab must be non-vanishing if there 
are to be gravitational waves in cosmology, thus refining the known covariant condition 
that the curls must be non-zero. Further investigation of the role of the distortions is 
clearly warranted. In particular, the covariant classification of cosmological solutions 
via E ab , H ab and their divergences and curls (see |J, [|J, f5|, 0, pTOf] , |i"7|), could be 
extended by considering also the distortions. 

Finally, we note that the electromagnetic analogue of these results follows from the 
covariant decomposition 

D b E a = D {b E a) + \D c E c h ab - |curl E c e ab c , 

and similarly for the magnetic field vector H a . This shows that the distortions of E a and 
H a are 'free' in the sense that they are not directly constrained by Maxwell's equations 
(see 0). Together with the curls, these distortions, as in the tensor case, are crucial for 
the existence of electromagnetic waves. 
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Appendix A. Covariant Ricci and Bianchi equations 

The Ricci identity for u a and the Bianchi identities (incorporating the field equations 
via the Ricci tensor), are covariantly split into propagation and constraint equations in 
PJ. For a perfect fluid source, these equations are in our notation the following: 



Ricci: 



+ |6 2 - D a ii a + \{p + 3p) = - a ab a ab + 2u a uj a + u a ii a , (Al) 

co {a) + |9w a + |curl?i a = a ab u b , (A2) 

&{ab) + fQ^ab — D(aMfo) + i?afe = ~ ^c{a^b)° 

- U( a uj b ) + ii( a Ub) , (A3) 

|D a 6 - D b a ab + curl oo a = 2e abc u b u c , (A4) 

V a tu a = ii a u a , (A5) 

H ab - curl a ab - B {a u b) = 2u {a u b) , (A6) 



Bianchi: 



p+(p + p)Q =0, (A7) 

(p + p)u a + D a p =0, (A8) 
E {ab) + 6£ a6 + \{p + p)a a6 - curl H ab = 3a c{a E b) c + uj c e cd{a E b) d 

-2u c e cd{a H b) d , (A9) 
#(a6) + QH ab + curl i? afe = 3<7 c ( a .£/ft} c + u c E cd ( a H b ) d 

+ 2u c e cd(a E b) d , (A10) 

D 6 £ o6 - |D a p = e afec A# cd - 3H ab uo b , (All) 

D 6 # a6 -{p + p)u a = - e abc a b d E cd + 3E ab u b . (A12) 



These equations generalise those given for the case of irrotational dust in || . In the 
case of small anisotropies and inhomogeneities, i.e. when the cosmology is almost FRW, 
the right hand sides of the above equations can be set to zero. The resulting linearised 
equations are the basis for covariant and gauge-invariant perturbation theory 

By projecting the Ricci identity for rank-2 tensors, we can generalise the identities 
in for the commutation of covariant time and spatial derivatives, and for the curl of 
the curl: 



Commutation of spatial derivatives (exact, non-linear): 

B [a D b] S cd = | (§0 2 - p) S [a ( c h b] V - 2S [a {c {E b f - \Qa b] d) ) 
+ 2h [a {c {E [b]e - §0a 6]e } S d ^ e - 2a [a {c a b]e S d ^ 
- 2S [a {c {uo b] io d ^ - u e uj e h b f + \Qe b] d \u e ) 
+ 2{a { J c e b]ef + a e[b e a] ^ f }^S d ^ 

+ V c {"b]Ue + &eb]efw f } S d)e - e abe uj e SW . (A13) 

Commutation of curl and time derivative (linearised): 

(curl S ab )' = curl S ab - §0curl S ab . (A14) 
Curl of curl (linearised): 

curl curl S ab = -D 2 S ab + §D (a D c ,S b)c + (p - ±0 2 ) ,S afe , (A15) 
where D 2 = D a D a is the covariant Laplacian. 
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